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DETAILED ACTION 
Drawings 

1 . The drawings are objected to because: 

On figure 3, blocks 102, 104, 106, 108, 114, 116, 118, and 120 must use 
descriptive language to identify each of such blocks. 

A proposed drawing correction or corrected drawings are required in reply to the 
Office action to avoid abandonment of the application. The objection to the drawings 
will not be held in abeyance. 

Claim Rejections - 35 USC §103 

2. The following is a quotation of 35 U.S.C. 103(a) which forms the basis for all 
obviousness rejections set forth in this Office action: 

(a) A patent may not be obtained though the invention is not identically disclosed or 
described as set forth in section 102 of this title, if the differences between the subject 
matter sought to be patented and the prior art are such that the subject matter as a whole 
would have been obvious at the time the invention was made to a person having ordinary 
skill in the art to which said subject matter pertains. Patentability shall not be negatived 
by the manner in which the invention was made. 

3. Claims 16, 18 are rejected under 35 U.S.C. 103(a) as being unpatentable over 
Jones et al. (US patent # 6,128,351), in view of Komatsu (US patent # 6,144,860). 

Regarding claim 16, Jones et al. et al. discloses: 

- an input data for receiving a known data signal (column 1 , lines 53 - 54); 
an input data for receiving an unknown data signal (column 1, lines 5 1 - 52); 
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- a data combiner to produce the composite output signal having discrete data 
signals, wherein each discrete data signal comprising at least a portion of 
unknown data (column 1, line 61) and at least a portion of the known data 
signal (column 1, line 68). 

While Komatsu discloses a signal power ratio input for receiving a power ratio 
signal indicating a ratio combining the unknown data signal and the known data signal 
(column 3, lines 18-32). 

Therefore, it would have been obvious to one of ordinary skill in the art at the 
time of the invention was made to combine Jones et al. and Komatsu to improve the 
communication system. 

Regarding claim 18, Jones et al. further inherently discloses that combiner is a 
frequency domain data combiner (column 1, lines 22-23. As the combiner for receiving 
OFDM modulated signals, it must be a frequency domain data combiner). 

4. Claims 20 - 24 are rejected under 35 U.S.C. 103(a) as being unpatentable over 
Jones et al. (US patent # 6,128,351), Komatsu (US patent # 6,144,860), further in view of 
figure 1 admitted as a prior art by applicant. 

Regarding claims 20 and 21, respectively, as followed by the limitations analyzed 
in claim 16, Jones et al. and Komatsu differ from the instant claimed invention that they 
do not show the step of comprising a serial to parallel converter having a plurality of 
outputs. 

However, figure 1 admitted as a prior art by applicant discloses the serial to 
parallel converter (102) to provide a plurality of outputs. 
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Therefore, it would have been obvious to one of ordinary skill in the art at the 
time of the invention was made to combine Jones et al, Komatsu, and figure 1 admitted 
as a prior art by applicant to improve the communication system. 

Regarding claim 22 and 23, respectively, as followed by the Umitations analyzed 
in claim 21, figure 1 admitted as a prior art by applicant further discloses an inverse 
discrete Fourier transform module (104). 

Regarding claim 24, as followed by the limitations analyzed in claim 23, figure 1 
admitted as a prior art by applicant further discloses a cycUc prefix adder (108). 

Allowable Subject Matter 

5. Claims 17 and 19 are objected to as being dependent upon a rejected base claim, 
but would be allowable if rewritten in independent form including all of the limitations of 
the base claim and any intervening claims. 

6. Claims 1-15 and 25 - 41 are allowed. The following is a statement of reasons 
for the indication of allowable subject matter: 

Regarding to the claimed inventions, the prior art of record fails to show or render 
obvious of a multi-carrier communication system, comprising: 

A transmission unit comprising: 

- A data input for receiving an unknown data signal; 

- A data input for receiving a know data signal; 

- A power ratio signal input for receiving a power ratio signal indicating a ratio 
for combining unknown data and known data signals; and 

- A data combiner coupled to the data input for combining the unknown and 
known data signals in accordance with the power ratio signal to produce a 
composite signal comprising discrete data signals, wherein each discrete data 
signal comprises at least a portion of the unknown data signal and at least a 



ApplicationyControl Number: 09/778,306 
Art Unit: 2631 



Page 5 



portion of the known data combined in accordance with the power ratio 
signal. It also has an output adapted to provide the composite output signal to 
a multi-carrier transmitter transmitted a transmit signal on a communication 
channel, wherein the transmit signal includes the composite signal; and 
A receiving unit comprising: 

- A multi-carrier receiver for receiving the transmitted signal on the 
communication channel, having an output for providing a corresponding 
composite signal, wherein the composite signal comprises corresponding 
discrete data signals and is shaped by at least one signal shaping characteristic 
of the communication channel; 

- A channel estimator having a known data input, an input coupled to receive 
the corresponding composite signal, an input coupled to receive the power 
ratio signal, and an input for receiving at least one estimate of the 
characteristic of the unknown data signal, estimating the at least one signal 
shaping characteristic of the communication channel from at least the 
corresponding composite signal, the at least portion of the known data signal 
and the at least one estimate of the unknown data signal. The channel 
estimator has also an output for providing at least one estimated 
communication channel characterizing signal; and 

- An equalizer coupled to receive the corresponding signal, the known data 
signal, the power ratio signal, and the at least one estimated communication 
channel characterizing signal, configures at least one of its signal shaping 
characteristics to compensate for the at least one signal shaping characteristic 
of the communication channel, shapes the corresponding composite signal 
accordingly. The equalizer also has an output for providing at least subsequent 
estimate of the unknown data signal. 

Conclusion 

7. The prior art made of record and not relied upon is considered pertinent to 
applicant's disclosure. 
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US Patent Documents: 

Barton et al. (US patent # 6,654,431 Bl) discloses a multi-carrier personal access 
communication system. 

Maloney et al. (US patent # 6,546,256 Bl) discloses robust, efficient, location- 
related measurement. 

Barton et al. (US patent # 6,449,246 Bl) discloses a multi-carrier personal access 
communication system. 

Jones et al. (US patent # 6,307,892 Bl) discloses a multi-carrier communication 
system and its corresponding method for peak power control. 

Humphrey et al. (US patent # 6,130,918) discloses a method and its 
corresponding apparatus for reducing the peak-to-average ration in a multicarrier 
communication system. 

Bauml et al. (US patent # 6,130,918) discloses a method and its corresponding 
apparatus for reducing the crest factor in digital transmission procedures. 

Bottomley et al. (US patent # 5,909,465) discloses a method and its corresponding 
apparatus for bi-directional demodulation of digitally modulated signals. 

Sakoda et al. (US patent # 5,907,583) discloses a transmitting/receiving apparatus 
and communicating method. 

Van Nee (US patent # 5,841,813) discloses a digital communication system using 
complementary codes and amplitude modulation. 

Other Publications: 

MioUsavljevic et al., "Fixed Point Algorithm for Bit Rate Optimal Equalization in 
Multicarrier Systems", IEEE International Conference on Acoustics, Speech, and Signal 
Processing, ICASSP 1999, vol. 5, 15 - 19 March 1999, pp. 2515 - 2518. 
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Wang et al., "Joint Channel Estimation and Equalization in Multicarrier 
Modulation System Using Cyclic Prefix", IEEE International Conference on Acoustics, 
Speech, and Signal Processing, ICASSP 1999, vol. 5, 15 - 19 March 1999, pp. 2733 - 
2736. 
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If attempts to reach the examiner by telephone are unsuccessful, the examiner's 
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ABSTRACT 

In this paper we present a fixed point algorithm for bit rate 
optimal time domain equahzation in multicarrier systems. 
The algorithm outperforms standard MMSE equalizers in 
capacity measure. Moreover, it adaptively estimates the 
optimal signal delay which was an exhaustive search prob- 
lem in MMSE equalizers. We briefly discuss the conver- 
gence properties and demonstrate the performance of the 
algorithm on a twisted pair copper wire loop with additive 
non-white noise. 

1. INTRODUCTION 

Multicarrer modulation (MCM) systems have become the 
technology of choice for high-speed data transmission over 
spectrally shaped channels. MCM transceivers divide the 
available frequency spectrum into a large number of paral- 
lel, independent and approximately flat subchannels, thus 
directly implementing the maximum hkelihood (ML) detec- 
tors on the receiver side. The preferred form of multicarrier 
modulation is the Discrete Multitone Transform (DMT), 
which uses the computationally effective Fast Fourier Trans- 
form to create independent channels. Combined with the 
optimization of bandwidth and bit distribution across us- 
able subchannels, this modulation scheme demonstrated 
the results comparable to the optimum eigendecomposition 
based MCM trainsceiver of [9, 3, 1]. 

The fundamental implement at ion al difficulty of the mul- 
ticarrier transmission is the intersymbol interference (ISI) 
which renders the subchannels mutually dependent. Inter- 
symbol interference is usually combated by addition of the 
cycUc prefix whose length is equal to the channel memory, 
so that the output sequence seems periodic to the channel. 
This quasi-periodicity makes the channel description matrix 
circulant. In case of additive white noise, the FFT basis 
vectors become the eigenvectors of the channel description 
matrix and its eigenvalues become the DFT coefficients of 
its first column. 

The addition of the cyclic prefix, although very efficient 
in decoupling the subchannels, does decrease the transmis- 
sion rate, especially for the channels with long finite im- 
pulse response (FIR) and small number of DMT carriers. 
The equalization in MCM systems is therefore divided into 
two parts: time domain equaUzation (TEQ) and frequency 
domain equalization (FEQ). TEQ focuses on the task of 
optimaly (in a certain sense) shortening of the impulse re- 



sponse of the channel to the length of desired cyclic prefix, 
while FEQ equalizes the frequency response of the channel. 
After the TEQ the cyclic prefix is discarded together with 
the intersymbol interference contained in it. 

Various approaches to the time domain equalization are 
possible. Due to ease of their implementation and exis- 
tence of adaptive algorithms, the minimum mean square 
error (MMSE) algorithms are most widely used. MMSE 
algorithms minimize the error between the target impulse 
response (TIR) and the equalized channel response. The 
length of the target impulse response is selected to be same 
as the desired length of the cyclic prefix. The minimization 
of the mean square error therefore results in minimization 
of the trailing ISI in the residual energy sense. The solution 
of the MMSE problem for DMT transceivers exists in the 
closed form [9, 3]. However, this structure does not guar- 
antee the optimal ity of the channel capacity, which is the 
approach we discuss in this paper. The channel capacity 
optimization algorithms have been discussed by Al-Dhahir 
and Cioffi [3, 1] in the form on nonlinear optimization prob- 
lem. Although it has been shown that this approach has a 
definite performance advantage to MMSE it has also been 
demonstrated that it is computationally very intensive and 
has very uncertain convergence properties. 

In this paper we analyze the problems behind the ca- 
pacity optimization in more detail and reveal why most 
common optimization techniques fail. However, the main 
contribution of this paper is introduction of a fixed point 
technique having following benefits: 

1. Guaranteed convergence to capacity optimizing val- 
ues of FIR TEQ equalizer regardless of the starting 
point, 

2. Adaptive determination of optimal signal delay, and 

3. No matrix inversion smd autocorrelation matrix esti- 
mation is required. 

The algorithm however, does have substantial memory re- 
quirements. 

2. PROBLEM ANALYSIS 

The block diagram of a DMT transmission system is shown 
in Figure 1. A iV-point FFT is used to divide the channel 
spectrum into N subchannels. For channels to be com- 
pletely independent and memory less, N has to be infinite. 
In practical applications, on the other hand, N is chosen 
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Figure 1: Block diagram of a DMT transceiver. 



of the probability of error, required margin and the coding 
gain. In order to treat all bins equally we set F,- = T. 
The overall number of bits can be expressed as: 



(3) 



Figure 2: Block diagram of time domain equalization. 



to be a power of 2. A common number used is 512, The 
finite number of subchannels results in the inter-block in- 
terference (IBI) which is mitigated by addition of the cyclic 
prefix (CP) which clears the channel memory thus mak- 
ing the successive block transmissions independent. Thus 
transformed signal is then low-pass filtered, converted to 
analog and sent through the channel. On the receiving side 
the procedure is reversed. 

As we indicated before, for highly dispersive channels, 
the addition of the cyclic prefix whose length is same as 
the length of the channel response severely decreases the 
data rate. The solution to this problem as suggested in 
[9] is to linearly equalize the channel impulse response to 
a shorter, target impulse response (TIR). Target impulse 
response is chosen to be closest in the mean square error to 
the channel-equalizer combination as in Figure 2. 

Framework we adopt is one of fractional equalization. 
The time domain equalizer is taken to be of length TV/ + 1 
and is denoted by vector w = [wo, wi, ... yWjVf]^ . Target 
impulse response is taken to be of length iVt + 1 and is 
denoted by b = [6o,6i, . . . ,6^;]^. 

For any fixed b, from the orthogonality principle we 
have that the filter coefficients of the feed-forward equalizer 
w can be written as: 



w*Ryy =b'R^ 



(1) 



where Ryy and R^y are the crosscorr elation matrices. This 
is the minimum mean square error (MMSE) solution that 
has been well documented in literature [9, 3, 1]. 

As it has been noted above, there is no guarantee that 
this solution will provide us with the maximum channel 
capacity In order to investigate the channel capacity prob- 
lem we compute the number of bits in each of iVfrequency 
bins of DMT and sum them to obtain the number of bits 
transmitted in one DMT symbol: 



(2) 



where V is the signal to noise ratio gap that characterizes 
the distance (in dB) between SNRi and the signal to noise 
ratio required to achieve the capacity and it is a function 



where we introduced SNR^.om.tric = V [n^lx ^ + ^] ' 
r. Typically, we assume that factors "-hi" and "-f" can 
be ignored and simplify the above expression to: 



SNR, 



'geometric ' 



(4) 



Since the value of capacity is directly related to the 
SNRgcometric, from hcrc on we consider only the optimiza- 
tion of this value. For the purpose of analysis we express 
the SNRg 

eometric lu tcrms of coefficicnts of vector b as: 



b^Mjb 
b'Nib 



where 




(5) 



(6) 



and, *i , ^J*""^ , and are the Fourier matrices such that 
= b^^f'^^b. \Wl{\\^ = w'-^r^'w and W[iYB[{\\ = 

Since quantity ^'^ ' ^ is recognized as the Rayleigh co- 
efficient, it follows that the vector b that maximizes the 
capacity for bin % is given as a solution of the generalized 
eigenvalue problem for the pencil (Mi,N|]. Unfortunately, 
this conclusion cannot be generalized for the product of 
Reyleigh coefficients for each of the frequency bins. It is 
not clear weather the closed form solution for this problem 
exists. 

Even the numerical optimization of the log-product of 
Reyleigh coefficients turns out to be cumbersome. First, a 
simple gradient optimization is computationally intensive 
since the gradient is given by: 



eometric = SNR, 



geometric 



b*M, 



b*Ni 



b*Nib b-Mib 



(7) 

and, therefore, requires quadratic computational effort. 

Moreover, due to the specific structure of matrices Mi, 
the numerator of the SNRgcometric has numerous zeros. For 
example, when iV^ = 4 (i.e. the cyclic prefix is of length 5) 
and for FFT length of 64, the numerator of SNRg^omrtric 
has 101 zero eigenvalues. Moreover, matrices can be 
poorly conditioned and may have the minimum eigenvalue 
close, or even equal to the small or zero eigenvalue of the 
numerator. OverEill, this structure contributes to the very 
hard optimization surface with many local minima and max- 
ima and possibly discontinuous gradients. 
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Figure 3: Dependence of capacity on the parameter 9 — 
arctan 62/61 for a simple channel. 



In order to better illustrate above described behavior 
we performed a simple exercise that can be easily tracked 
analytically. We selected the order of filter b to be 2, and 
used a simple channel: [1,0.9,0.9). We computed the ca- 
pacity for all values of the angle B that defines the unit norm 
vector b as 0 = tan" '{62/61). The parameter Q describes 
the whole space of vectors b, so the optimization surface 
is actually a line. The dependence of the capacity on the 
parameter Q is shown in Figure 3. 

It is obvious that for a more complex surfaces standard 
optimization algorithms suffer from locality problems and 
cannot guarantee the convergence to the optimal capacity. 

3. FIXED POINT HYBRID EVOLUTIONARY 
ALGORITHM 

3.1. Description of the Algorithm 

We propose an alternative way of DMT TEQ filter de- 
sign well suited for the fixed point DSP processor appli- 
cations. The algorithm is a hybrid of the Simulated An- 
nealing and Genetic Algorithms and by construction bor- 
rows the global convergence from Simulated Annealing and 
convergence speed from Genetic Algorithm. 

We initialize the algorithm with a set (population) of TV 
column vectors r — [b',w\ A]' and the control parameter 
T. A here stands for the received data delay required for 
the best symbol alignment. 

For each of the vectors we compute the capacity and 
rank and sort the vectors with respect to their capacity. At 
this point it is important to indicate that the capacity com- 
putation is not hard, nor computationally intensive since it 
can be performed via efficient hardware implementation of 
FFT. 

Prom the sorted population of vectors we choose pairs 
in a top-to-bottom fashion. These pairs become parents 
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Figure 4: Flowchart of the hybrid filter design algorithm. 



to new members of the population. Two new vectors are 
created from each set of parents by exchange of random 
number of bits whose position is selected at random as well. 
These new vectors are called children. The Boltzmaan trial 
is held between the parents and children in following sense: 

1. If the capacity of parent vector is smaller than the 
capacity of at least one of the children vectors, the 
better performing child vector takes the place of the 
parent. 

2. If the capacity of parent vector is greater than the 
capacity of both child vectors, child vectors take the 
place of the parent with probability: 

p = e (8) 

We then proceed to mutate the population. This is 
performed by randomly selecting the pre-determined ratio 
of the population, and altering the random bits in the vector 
representation. Finally as the last step in the algorithm, we 
reduce the control parameter T in a logarithmic fashion. 
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3.2. Convergence Issues 

Prom the description of the previous subsection it is ob- 
vious that the algorithm has strong similarities with the 
Simulated Annealing since it uses the probabilistic tran- 
sition rule governed by the exponential of the difference 
between performances of the two candidates. The devia- 
tion frotn standard Simulated Annealing is twofold: there 
is 7V/2 transitions happening in every iteration (generation) 
and neighborhood generation is dependent on the multiple 
members of the population. 

For the proof of convergence we use idea of [10] to con- 
catenate all vectors r in one large super- vector We define 
the performance criteria of this supervector as the sum of 
performances of individual vectors r it contains. The algo- 
rithm then becomes the Simulated Annealing algorithm of 
[11] over the extended space of supervectors f and asymp- 
totic convergence of supervectors readily follows. This, in 
turn, implies convergence of the population to the singular 
distribution at the maximum capacity point. 

4. EXPERIMENTS 

In order to demonstrate the effectiveness of the proposed al- 
gorithm we have implemented the algorithm for the CSA6 
9kft - 26ga copper wire test line. The transmitted signal 
power was -40dBm which is standard for the commercial 
DSL services. Additive noise was modeled as ETSIB. The 
algorithm was initialized with initial population of 5000 vec- 
tors, and run for 100 iterations (generations). The muta- 
tion factor was constant and equal to 20% of the population 
vectors and 20% of the bits in the mutated vectors. The 
mating was performed by exchanging 40 % of randomly se- 
lected bits from both parents and each parent was paired 
in the Boltzmann tournament with only the child closer in 
the Euclidean distance to reduce the computational load. 

The results of the experiment are shown in the Figure 
4. The proposed algorithm clearly outperformed the MMSE 
as early as generation 7 and converged at generation 40 to 
a value close to the maximum available capacity. 
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ABSTRACT 

Multicarrier modulation (MCM) is a promising technique for high 
rate data transmission. A one-tap equalizer is an essential part 
of the MCM system and the channel estimation is needed to get 
the coefficient of the equalizer. Lack of correct channel estimation 
may cause significant performance degradation. We propose to use 
the cyclic prefix to estimate the chaimel for the MCM system. We 
found that the cyclic prefix originally used solely to guarantee the 
optimahty of modulation using discrete Fourier transform (DFT) 
can be viewed as a source of channel information. Based on this 
observation, we propose a joint channel estimation and equaliza- 
tion algorithm using the cyclic prefix. Our simulations show that 
the algorithm can adaptively track the variation of a moderately 
time varying channel and has about l-2dB gain over the system 
using the charmel estimation obtained by the conventional training 
schemes. 



1. INTRODUCTION 

Multicarrier modulation (MCM) is now considered an effective 
technique for high rate data communications in both wire and wire- 
less environments. The principle of MCM is dividing the transnut 
data into several parallel low bit rate data streams, and using these 
data streams to modulate several carriers, which in frequency do- 
main is equivalent to partition the entire channel into several par- 
allel subchannels. 

MCM provides an optimal way for channel capacity usage by 
adjusting the bit rate and transmit power according to the SNR 
of subchannels. MCM also has a relative longer symbol duration 
since it is a block oriented technique. The long symbol duration 
produces greater immunity to impulse noise and fast fading. Be- 
cause of these advantages, MCM is considered a promising ap- 
proach in digital subscriber line (xDSL), digital video/audio broad- 
casting, and wireless communications. 

In MCM system, usually a one-tap equalizer is needed for 
each subchannel to get the estimations of transmitted data. The 
channel information is essential to the coefficients of the equal- 
izers. Some techniques, such as differential PSK modulation, are 
used to eliminate the need for channel estimation and equalization. 
However, differential demodulation causes 3-4dB Signal to Noise 
Ratio(SNR) loss compared with coherent demodulation if channel 
information is known. Moreover, channel information is also very 
important for the bit and power allocation. 

In appUcaiions such as xDSL. some training processes are per- 
formed to estimate the channel before the communication is set 



up. Then, this channe' estimate is used through the entire com- 
munication. If the channel changes, retraining is required to track 
the variation. Recently some research has been done on channel 
estimation and tracking in wireless communications. A minimum 
mean square error estimation algorithm is proposed in [5]. 

In ttiis paper wc propose a new channel estimation scheme that 
can track the change of the charmel parameters without retrain- 
ing. Usually a cyclic prefix or a guarding period is added between 
two symbols in MCM system in order to reduce the intersymbol 
interference (ISI). We propose to use the cyclic prefix, which is 
normally discarded, for channel estimation and equalization. We 
observed that the prefix actually provides a constantly sent training 
sequence if accurate transmit signal can be recovered by the con- 
ventional MCM systems. A joint channel estimation and equal- 
ization algorithm using the cyclic prefix is proposed based on this 
observation. Simulations were performed under the asymmetric 
digital subscriber line{ADSL) environment to show the effective- 
ness of the algorithm. 

Z MCM SYSTEM USING CYCLIC PREFIX 

MCM partitions a spectrally shaped channel into a number of par- 
allel and subchannels by modulating a set of orthonormal basis 
functions. Most of the MCM systems choose the inverse discrete 
Fourier transform (IDFT) as the orthonormal basis. Fig. 1 shows a 
MCM system using IDFT as modulation scheme. 

Input data are first buffered into blocks which are used to 
form the symbols transmitted in channel. Each block of data is 
then divided into m/2 bit streams in a nwnner determined dur- 
ing system initialization and mapped to some complex subsym- 
bols to form the input of a m-point IDFT which is represented as 
Xfe = [Xo.fc Xi^k ' ■ ■ A'm-i,*]^, where JCi.jt is the tth input of 
IDFT. The modulation is then performed by m-point IDFT and the 
resuU is x* = [xo,k ^i,* ■ • * aCm-i.fc]^. 

The channel is usually modeled as a FIR filter with length v -\- 
1. The impulse response of the channel is h = [fto, /ii , • * * , f^v]^- 
To reduce the ISI caused by the channel memory, a cyclic prefix 
x^^* = • • * a:_i,it]^, which consists the last v samples of 

xjt, i.e.. 2;_i,jfc = im-i,*, t = 1, ■ - • , V, is appended in front of 
x;fc before transmission. 

At the receiver, the prefix part y[^^ = [v-v.Jt * * - y-i.fc]^ is 
discarded, only = [j/o.jk yi,k ■ • • j/m-i,*]^ is used for demod- 
ulation. The demodulation is performed by the DFT operation and 
the result is Yjt = [Ko.Jt Vi.k • ■ • Ym-ukf. 

It can be proved that the above modulation scheme is optimal 
due to the use of cyclic prefix in the sense that the mutual infor- 
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Figure 1: MCM System with Cyclic Prefix and Adaptive Channel 
Estimation 



mation is maximized. As m goes large, the subchannels can be 
viewed as independent with each other, i.e.. 



V 



(1) 



are samples of m point DFT of h. Ni^k are samples of m point 
DFT of the channel noise. 

Assuming Ni^k are independent with each other, the best esti- 
mation oiXi^k from yi,jt is achieved by applying a one-tap equal- 
izer tWi to Yi^k, i.e., 

Xi^k = Yi^k • ii^i- (2) 
The optimal coefficient for the one-tap equalizer is: 



where Vi is the transmitted power of Xi.jt and Ai = E[||yVi,jfc||^]. 

Then, Xi^k is the hard decision result of X^.jtt i.e., Xi^k = 
q{Xi^k)> where g(-) is some kind of quantization function. 

3. THE EXISTING TRAINING METHOD 

The channel is modeled as the FIR filter stated before. When a 
training sequence xt is sent to a channel, the output of the channel 
is: 

where njk is uncorrelated random noise. Suppose hk is the esti- 
mation of hk, then the estimated output yjt ~ xk * hk. The best 
hk is chosen to minimize the power of the error ejt — yk — yk- 
This is the familiar quadratic form of minimizing the mean square 
error problem. There are many well-known methods to solve this 
problem, such as least squares (LS) method which is also used in 
our algorithm. 



For MCM systems, we need to estimate the channel parame- 
ters in frequency domain. Instead of the lime domain algorithm, 
an equivalent frequency domain deterministic least squares (DSL) 
channel identification algorithm can be used. In this algorithm, 
a training block with length m is sent periodically, and then the 
channel outputs are collected and averaged to reduce the influence 
of channel noise. The DFT of channel response is obtained by 
performing element by element division between the DFT of the 
averaged channel output and the input training sequence. Several 
different training blocks with guarding band can be used in order 
to further average out any non-linear effects. The final estimation 
is obtained by averaging the results of all these training blocks. 

4. THE PROPOSED JOINT CHANNEL ESTIMATION 
AND EQUALIZATION ALGORITHM 

4.1. Observation on Cyclic Prefix 

The training algorithms above is designed for the time-invariant 
system, which means new training process must be performed if 
the channel varies. However, in the MCM system using cyclic 
prefix, we can view the cyclic prefix as a training sequence and 
use it to track the variation of the chaimel. 

Let's consider the prefix part yj^^^ which is originally dis- 
carded. The relationship between y[^^ and the transmit signal is 



where 



A*: 



(4) 



The lower triangle part of matrix is composed by 
while the upper triangle part is composed by the last v~l samples 
of Xfc. However, this last t; - 1 samples are also the elements of the 
prefix x^{_\. So if all the prefix parts concatenate together as a pair 
of sequences x*^^ = {• • • x_„,fc_i • ■ x-i^fc-ix-v.* ■ • -x-i.t • • ■} 
and y^^^ - {- - y-v,k-i • • ■ y-i,k-\ y-v.jt ■ - - y-i.* • * the re- 
lationship between these two satisfies 



Vk 



if) _ ^if) 



(5) 



If we can get accurate estimations of transmitted prefix by the 
conventional MCM method, i.e., we know x^-^^ then (5) shows 
that y^^^ and x^^^ form a pair of training sequences that can be 
used to estimate the channel. 

One problem here is that we can only get the estimate of this 
training sequence. This estimate forms a feedback loop for the 
charmel estimation. The error incurred by the inaccurate estima- 
tion may propagate. In order to reduce the probability of error 
propagation the samples after hard decision, Xjt, are used to esti- 
mate the transmit prefix. Since the prefix is a time domain signal 
while Xfc are in frequency domain, an IDFT is performed to get 
the time domain estimation of xjt. 
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4.2. Least Square Method to Estimate h 

Several methods have been tried to solve h from the training se- 
quence formed by the cyclic prefix. The following method has the 
best performance in simulation. 

The idea of this method is trying to use LS method directly to 
solve (4). However, it is observed that Afc is an under-determined 
matrix. In order to reduce the efifect of random noise, we expand 
(4) to form the following equation: 



(6) 



After arranging data to the above form, the LS solution is given 
by 

whercAA^,L(A) = [Afc_;v^- ■ • Ajt • • • Ajt+L]^ and y^.i(Jb) = 
]yi-s • • ■ y*^* * • • xI+l] • A)v^^{fc) is the pseudo inverse of 
A^,l(^) that can be obtained by performing singular value de- 
composition (SVD) on A;v.l(A). which is an (JV + L + l)v x 
{v + 1) matrix. Usually N,L>Q, so the rank of A^f,L(A:) is con- 
strained by the length of the prefix, i.e., Rank(AN,L{Jt)) < u -f- L 
The roles of N is some kind of similar to forgetting factor. L 
is used to guarantee the amount of data is enough to get an accurate 
estimatioa 

43. Joint Channel EstimatioD and Equalization Algorithm 

Based on the discussion in section 4. 1 and 4.2, we summarize the 
channel estimation and equalization algorithm as the following: 

Input: received prefix part y^f^ and demodulated signal Y it- 
Known parameters: transmitted power Ti and noise power A;. 
Selecting parameters: N and L. 

Initialization:^ A = 0, an initial training is used to get the estima- 
tion of h(0). 
Computation: k = 1, 2, 3, • • • 



rfni(k~iy 



Xi,k = Yi,kWi(k-l) 



If A: = nL, where n is an integer, use Xi^k calculated above 
to form the matrix AN,z,(Jt - L). then. 

Hk)^X^^^dk~L)y^,L(k-~L); 

otherwise. h(A;) - h(Jt - 1). 



Here, what we present is actually a block recursive algorithm 
and the channel estimation is refreshed every L symbols. The sym- 
bol by symbol recursion is just the special case as L = 1. The rea- 
son for such a scheme is that this algorithm is a feedback scheme 
which combines channel estimation and equalization together. It 
requires more most recently data to keep on with the channel vari- 
ation. Our simulation shows that both JV and X should be chosen 
carefully to get the best performance, usually L > 1. 

5. SIMULATION RESULTS 

In our simulation, the transmit power of all the used subchannels is 
set to equal and fixed to 1. QAM signal is used in each subchannel. 
At first, some target error probability is preset. Then the bit is 
allocated by the following error probability constraint 

P,<4Q(*P) 

where di is the minimum distance between the signal points in 
QAM constellation of the ith subchannel. 
Initially the channel transfer function is 



m(D) = 



0.1 -|-0.8£>^ 



1 - 1.5D + 0.54J92 

The bit allocation is done according to this transfer function and 
will keep unchanged during the simulation. After some time, the 
channel transfer function will change to some H{D), Two differ- 
ent transfer function are used for H{D). 



Hl(D) ^ 



H2{D) = 



0.1 -h O.SD^ 
1 - 1.5D 0.54Z?2 

0.1 +0.8D^ 



1 - 1.4D -I- 0.5Z?2 * 
Length of FFT is chosen as m = 512. White noise is used in 
order to simplify simulation, i.e.. A; = A. 

The averaged mean square error (MSE) per subchannel is de- 
fined as 



err, 



1^1 

where em = \\Xi - Xi\\^ is the MSE of the ith subchannel and 
C/ is the set of all die used subchannels. ]U\ is the number of all 
the used subchannels. 

In our algorithm, first 2 frames of data are sent as pure train- 
ing sequence to get the initial channel estimation. After that, the 
real data arc sent and the joint channel estimation and equalization 
algorithm is used to track the variation of the channel. Fig. 2-4 
show the results of the simulation. The solid lines in Fig. 2 and 4 
show the result of the joint channel esdmation and equalization 
algorithm. The DSL channel identification algorithm is also per- 
formed for comparison. The dashed lines show the result of using 
the channel estimation obtained by DSL for HO{D) without re- 
training when channel changes, while the dash-dot lines show diat 
of using the channel estimation obtained by DSL for H{D) when 
channel changes. However, the training processes are not repre- 
sented in the following simulation results, since the block lengths 
of the training sequence and data transmission are different. It 
should be noted that only the channel estimations obtained by DSL 
are used in the following simulation and extra training sequences 
are needed in order to get those estimations. 
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(a) (b) 
Figure 2: Average MSE per Subchannel (A = 0.01, v = 128) 



Figure 3: MSE of the 88lh Subchannel (Pe = 10"'' , v = 128) 



In Fig. 2» the length of the prefix is 128. N and L are 2 and 
5 respectively. The average MSEs per subchannel are plotted. 
The channel changes from HO{D) to H1(D) in Fig. 2(a) while 
it changes to H2{D) in Fig. 2(b). We can see that the algorithm 
converges when the channel varies. However, it converges faster 
in (a), in about 10 symbols, than in (b), in about 100 symbols. If 
we consider that the channel change is more dramatic in (b) than 
in (a), the result is satisfiable. Moreover, the algorithm can not 
only track the channel variation but also achieve about IdB gain 
over the DSL method. In Fig. 2(a) the MSEs are plotted for both 
the target = 10""^ and = 10"^ In both cases, all 256 
subchannels are used and the SNRs of each subchannel are identi- 
cal. The only difference between these two cases is the minimum 
distance between the signal points, which means the errors of the 
estimation for the transmit prefix are different. The results shows 
that the adaptive algorithm is robust enough to such an estimation 
error. The MSE of f*e = 10"^ is only a slightly larger than that 
of Fe = 10~^. 

In Fig. 3, the MSE of the 88th subchannel are plotted for noise 
power A = 0.01 and 0.1 respectively, which means the SNR of 
this subchannel as A = 0.01 is lOdB higher than that as A = 0.1. 
This difference is compatible with the MSE difference in Fig. 3. 
As SNR goes down, the MSE is mainly caused by noise and the 
degradation brought by the inaccurate channel estimation becomes 
smaller. 

Fig, 4 shows the result with much shorter length of prefix, v — 
64. iV and L are chosen as 4 and 7. The other conditions of 
Fig. 4 are the same as those in Fig. 2 and the result is also similar. 
The performance gain over the DSL reaches 2dB which is even 
larger compared to the case with longer prefix length. However, 
the algorithm converges slower, in more than 100 symbols. This 
is because we need to collect more symbols to get an accurate 
estimation, i.e. since v is small, N and L must be large to maintain 
the dimension of AN,L(fc) to some level If and L axe too 
small, the channel estimation error becomes large enough to cause 
the error propagation through the equalization and the algorithm 
can not converge. As the prefix length goes even smaller, the error 
brought by the channel memory becomes too large that will also 
make the adaptive algorithm collapses. The smallest number we 
tried in our simulation that can still make the algorithm converge 
is 32. 




Figure 4: Average MSE per Subchannel (A — 0.01, v = 64) 



6. CONCLUSION 

We have presented a joint channel estimation and equalization al- 
gorithm using the cyclic prefix in MCM system. This algorithm 
can adaptively track variation of a moderately time varying chan- 
nel without additional training. Moreover, it also can give about 
l-2dB improvement over the conventional training schemes. 
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